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The influence of on-site (Hubbard) electron-electron interaction on disorder-induced localization 
is studied in order to clarify the role of electronic spin. The motivation is based on the recent 
experimental indications of a "metal-insulator" transition in two dimensional systems. We use both 
analytical and numerical techniques, addressing the limit of weak short-range interaction. The 
analytical calculation is based on Random Matrix Theory (RMT). It is found that although RMT 
gives a qualitative explanation of the numerical results, it is quantitatively incorrect. This is due to 
an exact cancellation of short range and long range correlations in RMT, which does not occur in 
the non-universal corrections to RMT. An estimate for these contributions is given. 
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I. INTRODUCTION 

The question considered in this paper is whether 
electron-electron interaction can reduce disorder-induced 
localization, thus enabling metallic behavior in two di- 
mensional disordered systems. 

The common view about the subject in the last 20 
years has been based on the well known scaling theory 
of localization^, according to which two dimensional sys- 
tems will always be localized (i.e., insulating), no matter 
how weak the disorder is. Although the original scal- 
ing theory did not take interactions into account, it was 
shown that weak interaction (i.e., high electron-density) 
does not affect its results^. On the other hand, in the 
limit of very strong interaction (i.e., very dilute sys- 
tems) it is known that the electron liquid freezes into 
a Wigner lattice, which is pinned by disorder and there- 
fore insulating 3 . All these results have lead to the opin- 
ion that the repulsion between electrons can only further 
decrease the conductance, so that all two dimensional 
systems will show insulating behavior, regardless of the 
strength of interaction between the electrons. 

A series of experiments performed in the last few years 
showed that even though in the limits of both very dense 
and very dilute systems we get the expected insulating 
behavior, for intermediate values of density (correspond- 
ing to r s between 4 and 40, where r s is the average 
inter-electron distance measured in the units of the Bohr 
radius) metallic- like temperature dependence is found 4 . 
The transition from an insulating behavior to a metallic 
one as the density decreases was entitled "Two Dimen- 
sional Metal-Insulator Transition" (2DMIT). An impor- 
tant feature of these systems is that an application of 
an in-plane magnetic field, (which cannot affect the elec- 
trons' orbital motion but can direct their spins) reduces 
the conductance in the metallic regime, until for high 
enough magnetic fields the conductance saturates as a 
function of the field, and the systems show the expected 
insulating behavior. This saturation field was estimated 
to be the field of full alignment of all the spins. 

These results arouse much interest and many ideas 



where suggested for their explanation. A debate started 
in the question of whether there is really a metallic behav- 
ior and a phase transition, probably caused by electron- 
electron interaction^; or the system is really insulat- 
ing, but the experimentally accessible temperatures are 
high enough to exhibit temperature dependent scatter- 
ing, thus causing the apparent metallic behavior 6 . 

Analytical^ and numerical^ calculations have shown 
that, as expected, for spinless electrons repulsion can 
only further localize the electrons, and does not lead 
to a metal-insulator transition. However, when taking 
spin into account, the situation is still unclear^. In a 
recent numerical exact-diagonalization study 8 , an An- 
derson model with both long range Coulomb interac- 
tion and short range Hubbard interaction was consid- 
ered. It was shown that the Coulomb interaction, exist- 
ing between any two electrons regardless of their spin, 
can only increase localization. On the other hand, not- 
too-strong Hubbard interaction were seen to cause de- 
localization (Strong Hubbard interaction will lead to a 
Mott-Hubbard insulator). Since this interaction exists 
only between electrons with opposite spins, its effect is 
decreased by an in-plane magnetic field, and disappears 
when all the spins are aligned. This dependence of lo- 
calization on interaction-strength and in-plane magnetic 
field thus mimics, at least qualitatively, the experimen- 
tally observed phenomena. Similar results were obtained 
recently using Quantum Monte-Carlo methods 9 . 

In this paper we wish to study further the weak short- 
range interaction regime, in which interaction-induced 
derealization was observed. We will first address the 
problem analytically, using a Random Matrix Theory 
(RMT) approach 10 , and then compare it to numerical 
simulations on an Anderson model. It will be shown that 
RMT can give only a qualitative but not a quantitative 
explanation for the numerical results, since RMT does 
not take into account non-universal correlations existing 
between wave functions in the diffusive regime. An esti- 
mate for the effect's order of magnitude and its depen- 
dence on the parameters of the system in the diffusive 
regime will be given. 
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II. ANALYTICAL RESULTS - RANDOM 
MATRIX THEORY 

We will consider an Anderson Hamiltonian with on-site 
Hubbard interaction : 



t ^2 & t;<r°-sr;<7 + U H ^2n S ;-[n s;l , (1) 



denote electron creation, anni- 



where a\. a: d sl - a and 
hilation and number operators, respectively, for a state 
on site s with spin projection a on some axis. The first 
term is a random on-site potential, where e s is chosen 
randomly from the range [-W/2,W/2]; the second is the 
hopping or kinetic term, where the sum is over nearest- 
neighbor sites s, st and t is an overlap integral; the third is 
the Hubbard term, the electrostatic interaction between 
two electrons in the same site (which must have opposite 
spins), whose strength is determined by the parameter 

To quantify localization, we will calculate the Inverse 
Participation Ratio (IPR), defined by P _1 = J2 S W^l • 
This quantity is of order 1 for localized states, and of 
order N^ 1 for dclocalized states, where N is the number 
of lattice sites. The IPR thus decreases when the single- 
particle wave function tp becomes less localized, and gives 
us an estimation for the changes in the conductance of 
the system. 

We assume here that without interaction the single 
electron energies and eigenvectors distributions for the 
ensemble of Anderson Hamiltonians are described by the 
corresponding distributions for an ensemble of Gaussian 
real symmetric matrices, i.e., the Gaussian Orthogonal 
Ensemble (GOE). This ensemble is defined by the well 
known distribution 10 : 
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where [3 = 1, A is a constant energy parameter, and fi(H) 
is a suitable measure. The eigenvectors are then a set of 
random orthogonal real normalized vectors. The average 
IPR without interaction for an electron in the n-th level 
with spin a is thus 11 
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where the superscript (0) denotes the state without in- 
teraction, and double angular brackets denote ensemble 
average. 

Now we add a weak Hubbard interaction, treating it 
in a self consistent way to first order in perturbation the- 
ory. Thus, the effect of spin-down electrons on the elec- 
trons with spin up will be the following effective potential 
(since only electrons with different spins interact, we have 
no exchange term): 
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According to the familiar first order perturbation the- 
ory, the first order change in the IPR of a spin- up electron 
in the n-th state due to its interaction with a spin-down 
electron in the m-th state is : 
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Since the wave functions can be chosen to be real due to 
time reversal symmetry, we omitted absolute value and 
complex conjugate notations in this and the following 
expressions. 

According to RMT, the eigenvectors distribution is in- 
dependent of the eigenvalues distribution, so we can sep- 
arate the averages of the numerator and denominator in 
the above expression. 

As for the average of the numerator, its value can be 
found in the literatureii* 1 ^, and the results are summa- 
rized in Table [I] We note that when s = s/ we have an 
average of even powers of wave functions at different sites, 
which is expected to be positive and vary as TV -4 , since 
we have eight wave function values in the expression, each 
of which goes as iV" 1 / 2 . On the other hand, when s =/= s/, 
it may appear at first glance that since we have an aver- 
age of odd powers of values of wave functions at different 
sites, which are uncorrelated, we should get zero. How- 
ever, we get in this case a nonzero negative value, going 
as N~ 5 . This result is due to correlations resulting from 
the orthogonality requirement on the eigenvectors. 

To understand this, we may note that squaring the 
orthogonality relation ^ s ipj(s)ipk{s) = for j ^ k and 
averaging, using the known resuhiii 
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we find that 
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for s ^ ,s/, i.e., if two different wave functions have the 
same sign on one site, from orthogonality they will tend 
to have opposite signs on another site and vice versa, 
hence the above nonzero negative average. 

As for the average value of the energy denominator in 
Eq. (J5J, in principle it might be possible to calculate its 
value using RMT. However, to estimate the leading order 
we will assume the spectrum is composed of equidistant 
levels, with mean level spacing A. 

Combining all those results together, we get, to the 
leading order in N, the following result for the change in 
the IPR of a spin-up electron in the n-th level due to its 
interaction with a spin-down electron in the m-th level: 
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TABLE I: Values of the average of the numerator in Eq. © for all the possible combinations of level numbers 1, m, n and sites 
s, s'. 
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where $(fc) is defined by: 

- 1 
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fe=i 

We observe that for m = n the correction is always neg- 
ative (for n in the lower half of the band), i.e., the in- 
teraction between two electrons in the same state tends 
to delocalize them, which is the only way to reduce their 
mutual interaction energy. For m ^ n the correction 
will usually be positive, i.e., electrons in different levels 
repulse each other, resulting in further localization. As 
can be expected, the former effect is larger than the lat- 
ter, due to the identity of the two interacting electrons' 
wave functions in the former case. However, the order N 
difference between the case m = n and the case m ^ n 
is caused by an exeat cancellation of the leading order 
dependence on N between the single short range (s = s/) 
term and all the N — 1 long range (s ^ si) terms in 
the latter case, which doesn't occur in the former. We 
will see below that this cancellation, together with the 
negative sign of the result for m ^ n, is correct only in 
RMT. 

Thus, if the lowest rcj_ levels are occupied by spin-down 
electrons, the total change in the IPR of a spin-up elec- 
tron in the n-th level is: 



The main features in the behavior of AP" 1 are as fol- 
lows : For n < n± the negative contribution of the spin- 
down electron at the same level n as the affected spin-up 
electron dominates the usually positive contribution of 
the other spin-down electrons. Therefore, AP" 1 is neg- 
ative, but decreases in absolute value when increases. 
For n > n±, there are spin-down electrons only in lev- 
els different from n, thus AP^ 1 is positive and increases 
when rij increases. At n = nj, there is a discontinuous 
jump of AP^ 1 . In both cases, since A ~ A^ 1 in real 
systems (although not in RMT), the effect is of order 
A^ -2 , if we keep the concentration of spin-down electrons 
constant. (We neglect here the logarithmic factor coming 
from the function $(n)). A plot of these formulas will be 
shown in the next section, where these expressions will 
be compared to numerical results. 

III. NUMERICAL RESULTS 

In this section we will examine results of numerical 
calculations and compare them to the analytical results 
discussed above. Two model Hamiltonians will be con- 
sidered : an RMT Hamiltonian and an Anderson Hamil- 
tonian. It will be shown that their results differ by an 
order of magnitude as well as in other characteristics. 
The theoretical predictions will be shown to agree with 
the former but not with the latter, and reasons for the 
discrepancy will be given. 

A. Random Matrix Hamiltonian 



AP„ 1 = We will first consider the change in the IPR for a true 

RMT Hamiltonian. Since we consider here only the weak 
interaction regime, instead of solving the exact many- 
24 Uh M _ ILLzi j (<j>( ?v — n) — <I>(n — 1)) body problem we simply diagonalize first the Hamilto- 

A y AT J nian without interaction, and then use the wave functions 

48 Uh to construct the effective potential, given in Eq. Q . This 

+ W~ ( ^ ~ n > ~ l )> ' n - n i'-> potential is then used to calculate the wave functions and 

the IPR with interaction. The applicability of this one 
loop Hartree-Fock approximation is justified by the fact 
24n| Uh that the change in P~ Y was found to be linear in Uh, as 

N4 A (*(JV - n) - $(n - 1)) required. 

48 Uh /jjj/ -|\ -|\\ The matrix size chosen was 408 x 408, and the elements 

j\[4 ^ \ \ n > \ n n l )) j n n l- were chosen according to the distribution law in Eq. J2J . 

(9) We have chosen A = O.lt, so that the mean level spacing 
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FIG. 1: Change in the IPR of a spin-up electron due to its interaction with spin-down electrons, according to RMT. The 
change is plotted as a function of the level number of the affected spin-up electron for different numbers of spin-down electrons: 
(a) n\ — 50; (b) nj, = 100; (c) nj = 150; (d) nj, = 200. In all the graphs the line indicates the theoretical formula, while the 
dots indicate the numerical results. The numerical results are averages over an ensemble of 5 x 10 4 realizations of 408 x 408 
RMT Hamiltonians. The estimated error approximately equals the width of the numerical results. Further parameters are 
given in the text. 



is A = 0.0196£, approximately equal to the spacing in the 
Anderson Hamiltonian, Eq. lJT]l. used in the next section 
(0.022t to 0.025t for W between 2.0t and 4.0t). The in- 
teraction strength Ujj was taken as l.Ot. The calculated 
quantities were averaged over an ensemble of 5 x 10 dif- 
ferent realizations. 

The numerical results for the change in the IPR vs. 
the level number of the affected spin-up electron due to 
its interaction with different numbers of spin-down elec- 
trons, are shown in Fig. together with the theoretical 
formula, Eq. 0. The theoretical formula was corrected, 
taking into account that the mean level spacing is not 
constant across the spectrum, but varies according to the 
semicircle lawifl, 

m=pw = 2^^ N - E2 > (10) 

where p(E) is the density of states. 

As can be seen, there is a good agreement between 
the numerical and the theoretical results. All the main 
features discussed at the end of the previous section can 
be clearly seen in the numerical data. 



B. Anderson Hamiltonian 

Now we will discuss the changes in the IPR for the 
Anderson Hamiltonian given in Eq. JIJ. The calculation 
was performed in the same method as was used for the 
random matrix Hamiltonian (i.e., one-loop Hartree-Fock 
approximation) . 

We have chosen a 17 x 24 lattice, corresponding to 
a 408 x 408 matrix. As for the RMT calculations, we 
took Uh = l-0t, while four values of disorder were used 
- W=2.0t, W=2.5t, W=3.0t and W=4.0t. The results 
were averaged over 10 4 realizations of disorder. 

First, in Fig. [3 the value of the IPR without interac- 
tion is shown for the four values of disorder, as well as the 
RMT value, Eq. We can see a difference here, as the 
Anderson model gives higher values (more localized) of 
the IPR than RMT. The effect is caused by non-universal 
(i.e., beyond RMT) corrections to the IPR and is more 
pronounced for higher disorder. The corrections for the 
IPR were calculated using supersymmetry techniques 1 ^, 
resulting in P^ 1 — Prmt ~ ff 1 ^^ 1 (where g is the di- 
mensionless conductance). We can also see, as expected, 
that the levels near the band edge have higher IPR, and 
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FIG. 2: The IPR for non-interacting electrons in the Ander- 
son model. The IPR is plotted as a function of the level num- 
ber. The lowest curve shows the RMT value, while the other 
ones are the Anderson model results for W=2.0t, W=2.5t, 
W=3.0t and W=4.0t, from lower to upper, respectively. The 
results are averages over an ensemble of 10 4 realizations of 
systems on a 17 x 24 sites lattice. The estimated error ap- 
proximately equals the width of the numerical results. Fur- 
ther parameters are given in the text. 



are thus more localized, than levels near the center of the 
band. 

Now we move to interaction effects in the Anderson 
model. The results are shown in Fig. |3 with the same 
occupation numbers as those chosen in the previous RMT 
calculations, for the four values of the disorder. As in 
RMT, the change in the IPR is negative for n < nj. and 
changes sharply (though not discontinuously) at n = n^. 
Nevertheless, it doesn't change its sign there. Moreover, 
the change in the IPR is larger by about an order of 
magnitude than the one found from RMT. Also, even in 
the range n < n;, it increases in absolute value, rather 
than decreases, when n± increases. All this is in contrast 
with Eq. © and the discussion following it. 

Another point is that the effect increases with disorder. 
This is seen by comparing AP^ 1 for the same level n but 
different values of W; or by observing that, for the same 
value of W, levels near the band edge, which are more 
localized, show larger AP,^ 1 . 

The reason for these differences is the above mentioned 
cancellation between long range and short range wave- 
function correlations in RMT. As has been seen in our 
RMT calculations (Table QJ , the average of wave func- 
tions product appearing in the numerator of Eq. is 
of order N~ 4 and positive when the two sites considered 
coincide, but are only of order TV -5 and negative when 
the sites are different. Since there are N — 1 terms of the 
latter type for each term of the former type, their total 
contributions are of the same order but their signs are op- 
posite. Due to the equallity of the numerical coefficients 
of the two types of terms when the interacting electrons 
are in different levels, they cancel out exactly to the lead- 
ing order in N , leaving behind a small negative term, of 



order A^ -5 . Therefore, in RMT interaction between elec- 
trons in different levels increases their localization, oppo- 
site to the situation for electrons in the same level. From 
this followed the decrease in the absolute value of AP,^ 1 
as n j increases in the range n < nj, , its positive value for 
n > ni, and the overall A^ -2 dependence of the effect for 
constant density of spin-up electrons. 

All this is correct when g is infinite. For finite g there 
exist non-universal corrections to the wave-function av- 
erages. Those corrections were not calculated before for 
the averages required here, but their behavior can be 
conjectured from known corrections for simpler averages 
(like those in Eqs. (0 J7J 14 ). We may expect them to 
have the same N dependence and sign as the RMT value, 
but to be smaller by a factor of g. The corrections for 
the short range (s — sf) terms and long range (s ^ sf) 
terms will not, in general, have equal numerical coeffi- 
cients, even when the interacting electrons are in differ- 
ent levels. Hence, after summation over sf we are left 
with an order g^N^ A contribution instead of the order 
A^ -5 contribution in RMT. For this reason, although the 
non-universal corrections are of order g~ x , for most of the 
averaged terms they are about N times larger, so they 
will determine both the magnitude and the sign of the 
interaction-induced change in the IPR. Since the correc- 
tions for s — sf will, in general, have a long range part, 
persisting for s ^ sf and having the same sign for neigh- 
boring sites (although for larger distances we may expect 
some oscillations), their sign will dominate the overall 
sign of the results. We will thus get a negative change in 
the IPR not only from interaction between electrons in 
the same level but also when the interacting electrons are 
in different levels. Hence, AP^ 1 will always be negative, 
as can be seen in the numerical results. 

Moreover, repeating the calculations with the non- 
universal correction to the averages of wave functions 
product, we can estimate the dependence of the effect 
on the system parameters. We expect the total change 
in the IPR of a spin-up electron due to its interaction 
with ni spin-down electrons to vary as 



This expression does not include a factor coming from the 
sum over energy denominators, which has only a weak de- 
pendence on N and nj, (logarithmic for equidistant levels, 
a weak power law for a non-constant density of states). 
Because wave functions corresponding to neighboring lev- 
els are more correlated than wave-functions correspond- 
ing to far away levels, there is also a factor, which changes 
sharply (though not discontinuously) when we pass from 
n < n l to n > nj,, as seen in the numerical results. Since 
A ~ JV — 1 in real systems (although not in RMT), the 
effect is of order g~ 1 N~ 1 , if we keep the concentration of 
spin-down electrons constant. This is in contrast to the 
N~ 2 dependence in RMT. Because N/g is much larger 
than unity in our numerical calculations, we can now un- 
derstand the order of magnitude difference between RMT 
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FIG. 3: Change in the IPR of a spin-up electron due to its interaction with spin-down electrons in the Anderson model. The 
change is plotted as a function of the level number of the affected spin-up electron for different numbers of spin-down electrons: 
(a) n i = 50; (b) n l = 100; (c) n l = 150; (d) n L = 200. In all the graphs the curves correspond to W=4.0t, W=3.0t, W=2.5t 
and W=2.0t, from lower to upper, respectively. The results are averages over an ensemble of 10 4 realizations of systems on 
a 17 x 24 sites lattice. The estimated error approximately equals the width of the numerical results. Further parameters are 
given in the text. 



and Anderson model results. Thus, all the features of the 
numerical data can be explained by taking non-universal 
corrections into account. 

As we have mentioned before, the non-universal part of 
the IPR without interaction, i.e., the difference between 
the value of the IPR without interaction in the Anderson 
model and its value in RMT, varies as g~ l N~ l . Accord- 
ing to the our estimate, the change in the IPR due to 
interaction in the Anderson model also goes as g^N^ 1 . 
Thus, their ratio, AP- 1 /(p- 1 - P^ T ), should be in- 
dependent of g, i.e. of the degree of disorder. It should 
also be independent of the number of lattice sites N if 
the densities of spin-up and spin-down electrons are kept 
constant. Thus, this ratio may be used to test our con- 
jecture for the parametric form of AP~ X . 

We first test the g indepence of the ratio AP^ 1 / (P -1 — 
Prmt) by plotting it in Fig.^for systems with identical 
lattice sizes (taken to be 17 x 24, as in the previous calcu- 
lations), but different values of disorder. We can clearly 
see that the differences between curves corresponding to 
different W values are much smaller than the correspond- 
ing differences in Fig. [3| The only exception is the value 
W=2.0 (the lowest curve), which shows a marked differ- 



ence from the other W values. This is probably due to 
the fact that for W=2.0 disorder is not high enough, so 
the electrons' motion is not fully diffusive, and ballistic 
boundary effects may be important. 



We now test N independence of the ratio 
AP- 1 /(P- 1 - P R l IT ) by plotting it in Fig. for 
systems with the same value of disorder (taken as 
W=4.0) but different lattice sizes - 8 x 13, 13 x 19, 
17 x 24. In all the cases the densities of spin-up 
and spin-down electrons are approximately equal (the 
horizontal axis is not the level number of the affected 
spin-up electron as before, but the filling i/, defined 
as the ratio of the number of spin-up electrons n and 
the total number of lattice sites N). We can clearly 
see that the different curves are almost identical. The 
only exception is the small 8 x 13 lattice, whose slighly 
different behavior can again be attributed to ballistic 
boundary effects. 
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FIG. 4: Ratio between the change in the IPR of a spin-up electron due to its interaction with spin-down electrons in the 
Anderson model and the non-universal part of the IPR without interaction. The ratio is plotted as a function of the level 
number of the affected spin-up electron for different numbers of spin-down electrons: (a) raj, = 50; (b) nj. = 100; (c) nj. = 150; 
(d) rij = 200. In all the graphs the curves correspond to W=2.0t, W=2.5t, W=3.0t and W=4.0t, from lower to upper, 
respectively. The results are averages over an ensemble of 10 4 realizations of systems on a 17 x 24 sites lattice. The estimated 
error approximately equals the width of the numerical results. Further parameters are given in the text. 



IV. CONCLUSIONS 

In conclusion, we have shown how a spin-dependent 
interaction can cause derealization, at least for weak 
short-range interaction. Localized electrons highly re- 
pulse each other, especially if they have the same orbital 
wave function and thus a different spin. This results in 
a tendency for interaction-induced derealization. The 
effect on an electron with a given orbital level and spin 
direction is stronger if the same orbital level is occupied 
by an electron with an opposite spin, and increases with 
the total number of electrons with opposite spin. The de- 
localization is thus reduced by an in-plane magnetic field. 
All this is in accordance, at least qualitatively, with re- 
cent experimental findings' 4 and numerical simulations^, 
regarding the in-plane magnetoresistance. 

We have also seen that the main difference in the influ- 
ence of the Hubbard interaction between realistic finite 



g systems and the RMT stems from exact cancellation 
of the leading order long range and short range terms in 
the former. Thus, while in RMT a state is correlated 
only to the same state with an opposite spin (except for 
weak anti-correlations with all other states), for finite g 
correlations between different states lead to a stronger 
repulsion between these states resulting in a stronger de- 
localization due to the on-site interactions. Nevertheless, 
the order of magnitude and parametric dependence of the 
IPR can be calculated using RMT, once the non-universal 
corrections are properly taken into account. 
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